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abstract 


A  study  of  turbulent  filci  condensation  was  first  made  based  on  the 
two-phase  model  utilizing  the  eddy  kinematic  viscosities  of  Deissler  and 
von  Karman,  It  has  been  shown  that  the  prediction  of  the  average  heat 
transfer  coefficient  is  still  higher  than  the  liquid  motal  data  of  Misra 
and  Bonilla.  Under  the  promise  that  an  upward  vapor  fl.ow  exists,  it  has  been 
shown  qualitatively  that  '.n  upward  vapor  flow  can  greatly  reduce  the  Iieat  transfer 
coefficient.  Ihe  observed  scatter  of  data  may  be  attributable  to  the  variations 
of  vapor  velocity. 

This  report  has  been  reviewed  and  is  approved* 
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Introduction 

♦ 

In  a  recent  paper  using  Nusselt's  model,  the  author  has  made  an' 
investigation  of  turbulent  film  condensation  by  utilizing  the  eddy  kinematic 
viscosity  expressions  of  Deissler  and  von  Kaniian*  A  few  important  conclu¬ 
sions,  which  are  relevant  to  liqv  >>  metals,  were  that  (i)  the  heat  transfer 
coefficients  (with  respect  to  Reynolds  number)  approach  a  limit  as  Prandtl 
number  becomes  small,  (ii)  tlxe  limit  is  not  too  much  different  from  Nusselt*s 
laminar  case,  and  (iii)  the  results  are  rather  insensitive  to  the  value  of 
the  ratio  of  eddy  thermal  diffusivity  to  eddy  kinematic  viscosity.  Conclusion 
(iii)  is  a  collorary  of  (i)  which  implies  that  the  magnitude  of  the  molecular 
thermal  conductivity  is  dominant  over  that  of  the  eddy  conductivity  in  the 
small  Prandtl  number  range.  Conclusion  (ii)  is  not  very  encouraging  from 
the  practical  standpoint  because  the  experimental  data  are  much  loxrer  than 
the  prediction  based  upon  Nusselt's  laminar  model.  Therefore,  it  was  men¬ 
tioned  in  that  paper  that  a  proper  description  of  the  problem  must  in¬ 
clude  the  following  modifications:  (a)  non-vanishing  interfacial  shear  stress, 
(b)  inertia  effects,  (c)  convective  heat  transfer  and  subcooling,  (d)  inter¬ 
facial  waves  and  rippling,  (e)  interfacial  thermal  resistance,  and  (f)  contact 
thermal  resistance. 

The  purpose  of  this  paper  is  to  investigate  the  turbulent  condensation 
by  including  the  first  three  modifications.  The  problem  can  be  formulated 
within  the  framework  of  a  two-phase  boundary  layer  treatment,  and  the  inter¬ 
facial  shear  stress  is  that  caused  by  a  quiescent  vapor  at  infinity.  This 
type  of  problem  has  already  been  worked  out  by  others  [1,  2,  3]  but  only  for 
the  laminar  case.  Therefore,  the  feature  of  the  present  analysis  is  to  inc- 
t^umbers  in  brackets  refer  to  the  References  at  the  end  of  the  paper. 
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lude  the  effect  of  turbulent  transports  In  a  two-phase  film  condensation 
model.  A  discussion  will  also  be  presented,  which  will  offer  an, explanation 
for  discrepancy  between  the  theoretical  predictions  and  the  experimental 
data. 


Formulation 

V  ' 

The  physical  model  for  film-wise  condensation  is  as  shown  in  Figure  1» 
where  the  x-axis  is  parallel  and  the  y-axis  is  perpendicular  to  the  plate. 
Suppose  a  semi-infinite  plate  is  immersed  in  a  quiescent,  saturated  vapor. 
If  the  plate  is  hold  at  a  lower  temperature  t„  than  that  of  the  vapor  tg» 
the  vapor  will  condense  and  a  film  of  condensate  will  flow  downward  due,  to 
gravity.  At  the  interface  of  the  condensate  and  vapor,  a  motion  will-  be 
induced  in  the  vapor  duo  to  a  shearing  action  of  the  condensate  and  thus  a 
vapor  boundary  layer  will  be  formed. 

Under  the  assumption  that  the  variation  of  physical  properties  can  be 
ignored  within  the  temperature  range  involved,  we  can  write  the  following.  , 
boundary  layer  equations: 
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Vapor: 

(4) 

(5) 

vrhere  i),  and  K.  denote  the  total  kinematic  viscosity  and  total  thermal  dif- 

\ 

fusivity.  ife  shall  impose  the  following  conditions: 

u  =  V  =  0  and  t  =  t  ,  at  y  =  0 

t  =  t  '  at  y  =  S 

s 

u  =  V  =  0  and  t  =  t  at  y  =  oo 

G 

Since  a  constant  temperature  will  be  assumed  in  the  vapor,  the  equation  of 
heat  eiicrg;^'  does  not  appear  in  that  layer. 

At  tiio  interface,  we  also  demand  the  continuity  of  velocity,  shear 
‘trof..:.,  and  mass  ’lux,  as: 


u(S)  =  u(S) 

Ou/ay)j  =  (iu/iy)g 

L 


As  a  useful  approximation,  it  has  been  shovm  ^1,  3]  .that  the  fulfillment  of 
the  above  interfacial  conditions  can  be  achieved  satisfaptorily  by  simply 
considering  the  following  relationship: 

^  ^  'I 

,6)  u(5)_  Udy  *  ■»t(— )5’0 

Therefore,  without  solving  u  explicitly.,,  this  allows  us  to  investigate  the 
tvio-phase  condensation  by  treating  only  tli  equa^tions  in  the  condensate  layer. 
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Method  of  Solution 

.  The  condensation  problem  can  be  characterized  bjy  introducing  an  in¬ 
tegral  relation  which  states  that  the  increase  in  the  condensate  flow  pate  , 
is  entirely  due  to  the  vapor  being  condensed: 

■  ■  V  S  ^  S  t 

(7)  —  (  u  dy  -  .0 

Cp  o  X  ■'  o  d  y  6 

Since  conditions  (6)  and  (7)  are  already  in  an  integral  form,  it  is  con¬ 
venient  to  convert  (2)  and  (3)  into  a  similar  form.  If  equations  (2)  and  . 
(3)  are  integrated  in  y  from  8  y  and  conditions  (6)  and  (7)  are  incor¬ 
porated,  we  obtain: 

(8)  2{  u(iu/&x)dy  ♦  u|  (du/dx)dy  ♦  u^(8)(d5/dx)  -  g(5-y)  -^^Ou/dy)  ■  0 

y  o 

15  S  3^ 

(0)  I  u(d0/dx)dy  -(  (l--8)(3u/3x)dy  -  (1-0) [  (du/dx)dy - 
'y  y  o 

8 

-(  X/Cp  At)  I  I  (3u/dx)dy  ♦  u(6)(d8/dx)  |  ♦  /r^(d0/dy)  -0 

o  ■ _ 

where  0  =  t  -  t^/At»  The  use  of  continuity  equation  (1)  has  been  made  in 
eliminating  v'e  in  the  above,  .  ' 

Following  the  spirit  of  a  phenomenological  approach,  the  total  trans¬ 
port  coefficients  are  assumed  to  be  the  sum  of  two  contributions  -  molecular 
and  eddy.  As  in  the  previous  work  C4],  we  shall  adopt  the  eddy  kinematic; 
viscosities  of  Deissler  and  von  Karman  and  let  the  ratio  of  eddy  thermal, 
diffusivity  to  ^dy  kinematic  viscosity  be  a  constant,  'Ihe  assumption  of  a 
constant  ratio  may  be  criticized  because  of  its  dependency  on  the  feclet  > 
number  and  the  flow  geometry.  However,  this  seinres  for  our  present  purpose 
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of  atte.apting  to  assess  a  gross  effect  of  turbulent  transports  -  more  will  , 
V)e  said  later, 

Near  the  plate,  we  have  Deissler’s  expression  for  eddy  kinematic  vis- • 
cosity  as: 

♦ 

(10)  u  y  (1  -  exp(-n^uy/7  ))  ,  0  *  y  ^  y* 

and  v;e  have  von  Karman's  expression,  which  is  valid  at  a  distance  from  the 
platen  as: 

(11)  ^2  “  ^  » y*  *  y 

\7herc  n,  K,  and  y*  are  empirical  constants# 

In  the  laminar  T?,  and  fC  are  constants}  therefore j  a  sxniilarity 

t  t 

transfonrdtio'i  r-oNaoes  ..a.uai/lons  (B)  and  {<))  into  a  one-dimensional  problem. 
■:lo.-.eve:,',  la  our  case,  such  a  simplification  does  not  exist  and,  in  general, 
if  and  /f .  will  introduce  further  non-linearity  into  the  system.  This  sug- 
gests  the  necessity  of  treating  the  nuraerical  solution  as  an  initial-value 
problem.  In  doing  so,  it  is  advantageous  to  make  a  transformation  such  that 
the  domain  of  u  and  0  becomes  rectangular  instead  of  a  wedge-shaped  boundary 
layer.  This  can  be  done  by  introducing  new  coordinates,  y/S  and  ^ =  x  , 
then  equations  (8)  and  (9)  become: 

(12)  2j\ou/d5)dn  ♦  U +  (5VS){(  <*n  1  “ 

-  g(i  -  n)  ° 


V/ 


(13)  (\(Ji/»s)d,i*fx(ii)/i>SMl  * 

♦  (8'/5){xl'*»  *1  ■'*  '*1 )  *  *  ® 

0  H 

,  «•_  /ti.  1 ')  _  ^  and  S’  =  dS/d^. 

where  X  -  U  -  ^  “  - -  * 

C  At  .  . 

P 


MiiTierical  Procedure  .  ■  , 

In  a  plane  of  t  andn,  the  set  of  equations  (12)  and  (13)  ,<!«>  Ix! 

eenverted  into  a  system  of  partial  differential  equations.  This  is  done  ty 
first  setting  up  a  tuo-d'imensional  mesh  in  A?  and  4^  (see  Flgur 
Since  0  ^  1  -  5  .  the  total  number  of  n-mesh,  be  assigned  ab 

initio.  If  the  integrals  in  (12)  and  (13)  are  approximated  V  the  trapa- 
soidal  rule,  each  equation  gives  rise  to  a  system  of  S  partial  differential 
equations  -  aU  for  the  total.  That  is,  equations  (12)  and  (13)  can  be 

vrritten  as;  ^2 

2  4n£<r%(i»»/S!)i  *  B„4nE<T(»0/n)l  *  (d8/dt)«‘V'«)|E  ‘ 

(14)  n  °  . 

o  ’ 

(15)  4,  E(r«i(M/»t)i  »  ‘1  Eff  »bM)i  t  1n*1 1'f  ♦ 

♦  (d5/d^)frV5){T„2ff  Bi  *  tfftl  Bi}  « 

where  r  is  l/2  for  the  first  and  last  terms  of  a  sum  and  1  for  all  others. 

For  (14)  there  are  IM  unknowns,  (W>t)i,-.  (»tl/»<)K  >«> 

(»a/>^)  vanishes.  On  the  other  hand,  since  (»T/»t)o  and  (>T/>{)||  vanis  , 
there  L  only  »-l  unknowns,  (>T/»t)i,-.  Ol/Hh-l-  “5)-  .Therefor,, 

we  have  a  consistent  system  of  2B  partial  differential  equations,  each  of 


which  is  of  the  following  form;  , 

(16)  ^  ^ 

where  f  denotes  Uj  Tj  or  6« 

With  oonditions  at  S  =  this  is  ah  initial-value  problem  of  tracing 
the  5 -evolution  of  f.  Since  equation  (16)  uas  originated  from  the  parabolic 
boundary  layer  equations,  it  can  be  considered  as  a  psuodo-diffusion  type 
equation  uith  a  variable,  non-linear  diffusion  coefficient.  An  implicit 
scheme  will  then  be  used  in  writing  (16)  in  a  finite  difference  form,  as! 


-  Vl  /*■>  ) 


From  the  standpoint  of  convergence  and  stability,  the  implicit -scheme  . 
ia  considered  to  be  superiot  to  the  explicit  one,  but  this  is  always  achlev^l 
at  the  expense  of  extra  labors.  However,  this  is  not  true  in  our  Case  be¬ 
cause  equations  (14)  and  ,(15)  must  be  treated  as  a  system  Of  equations,pre- 
ferably  linear.  At  any  rate,  no  additional  work  has  been  incurred.  In  the 
interest  of  simplicity, '  equation  (17)  was  linearized  by  taking  rather, 
than  an  Average,  e.g.  (I*  +  )/2,.  ln''ol''«»  interatlons.  The  Improve 

ment  in  the  final  results  dW  not  justify  the,  added  oonpllcaUons.  ’  :  .y,; 


f ) 

o 


The  final  finite  difference  equations  based  on  an  implicit  scheme  are: 

k  k+1  N-1  ic  jj+i  ic+i  k  ,  ,  k+1  k+1  ^ 

(10)  U„  Un  +  2  S  Ui  Ui  +  %  Un  ♦(Un  /2){2  C  J* 

n+1  1 

\  t  v2  /  It  x2  /  It  -2  k  .  k  k.  % 

+  8  (1/2  8  )  {  (Un  )  +  2  E  (Ui  r  +  (%  )  +  Un  (2  2:  Ui  ♦  Un)}  ♦ 

n+1  1 

*(1/5“  )2(-‘V2'‘qb(Vt)^  (Ci  -  Ci  )  -  8  (1  -  n-MXMAn)  ♦ 

*  (3/2)  {  (oj;  )='  *  2  (Hi  )*  *  (h5  )‘]  *  (s'e  uJ  *  uJ  ) 

■  n+1  1 


h“  *  2  ^  ^  nr  *  ^  2  XiT"  *  * 

n+1  n+1 

*r“{2  2  *  ujs }  *  (i/5'')'(AfA,b(A>^  (Ci  -  )  * 

*  s''*Mi/6'‘){^  (2  p  hJ  *1^  )*r^oJ  • 

1  n***! 

•  <  ^  ♦  2  £  tJ  *  h“  t|S  .  2T^(2  "n  oj  *  ^  }* 

n+JL  1 

««  • 


*  21tJ  *  2  £  H^  ♦  oJ  J 


For  our  problem,  any  attempt  to  assess  the  a  priori  convergence  by  means  of 
an  amplification  matrix  f7]  does  not  seen  very  practical.  Therefore,  we  shall 
resort  to  a  trial-and-error  approach  in  that  sufficiently  small,  yet  not 
impractical,  and  n»j  will  be'  chosen  by  testing  the  numerical  results.  This 
is  also  supported  by  Lax's  theorem  which  states  that  the  stability  is  the 
necessary  and  sufficient  condition  for  convergence.  The  overall  error  bound 
for  the  present  method  can  be  estimated  as  0(A^)  +  0(n’'j).  .  ■ 

Equations  (1<3)  and  (19)  form  a  system,  of  2IJ  linear  equations  in  2N 
unknoxms,  Tj,***,Tj^  8,  at  the  forward  step  of^  ;  and  the  aug¬ 

mented  matrix  can  be  determined  from  the  information  available  at  the  back-  ' 
ward  step.  > 


Numerical  Results  • 

The  system  of  (18)  and  (19)  was  first  studied  for  the  laminar  case 

by  setting  and  =/f.  The  forward  marching  numerical  solution 

was  commenced  at  5oC“  cm)with  appropriately  chosen  initial  conditions.* 

It  is  not  the  purpose  of  this, paper  to  re-coinpute  the  laminar  case;  however, 

it  provides  an  excellent  check  on  the  correctness  and  accuracy  of  the  present 

numerical  method.  This  was  done  by  exploring  the  simil^irity  of  boundary 

layer  equations;  that  is,  a  particular  J-dependenqy  exists  so  that  certain 

quantities  can  be  made  invariant  inj  .  For  instance,  U/(g/»^),  T,  and 

S/5  (where  S  =  C^(C  At/x)x/P  (g/i»^)  J  is  Nusselt’s  laminar  film  .  ■ 
o  .  •  p  -  P  .  r 

thickness)  must  be  independent  of J ,  Indeed,  for  the  worst  case,  the  con¬ 
stancy  of  the  first  4-digits  for  U/(g/i;),  5-  for  T,  and  3-  for  ^/^q»  were 
observed.  In  addition,  the  laminar  results  will  serve  as  the  initial  con-  ■. 
ditions  for  the  turbulent  case,  .  . 

In  this  work,  we  shall  restrict  ourselves  to  computations  for  two 
Prandtl  numbers,  0,008  and  0.003,  which  are  typical  of  liquid  metals.  From  ■ 
the  numerical  standpoint,  the  use  of  (10)  and  (il)  presents  some  difficulties. 
First,  von  Karman's  expression  requires  the  first  and  second  derivatives  of  . 
.U  and  the  numerical  differentiation  alvjays  amplifies  the  roughness  of  data. 

As  must  vanish  once  due  to  the  upward  '(negative)  interfacial  drag,  its 

magnitude  also  varies  in  a  wide  range.  Another  distracting  feature  was' that 

^  First  the  numerical  solution  was  started  with  the  initial  conditions 
corresponding  to  Nusselt’s  resultft  Ely  similarity  as  the  solution 
progresses  in  J  ,  U/(g/i>'^),  T,  and  S/8q  approach  limiting  values  which 
are  .  (Jiff erent  from  the  initial  conditions.  The  limiting  values 

will  then  be  used  as  the  initial  'conditions  for  the  final  computations* 
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the  stability  seemed  to  be  sensitive  to  the  irregularities  resulting  from 
numerical  differentiations.  By  trials,  it  was  found  that  the  consistent 
and  conservative  derivatives  can  be  obtained  by  smoothing  the  data  with  a 
second  degree  least-square  fit.  Second,  the  values  of  eddy  kinematic 
viscosities  of  (lO)  and  (ll)  are  not  usually  identical  at  y  ;  therefore, 
there  is  a  jump  at  that  point.  This  is  inherent  in  the  present  numerical 
method  because  such  a  difficulty  can  be  mended  rigorously  in  Nusselt's 
turbulent  case  [4],  Therefore,  we  have  Introduced  the  following  artifice 
of  smoothly  connecting  (10)  and  (11)  J  ,  . 

(20)  »*  -  »*  «tp(-0.695(y/y»)b  *  *2  (1  -  .xp(-0.69S(y/3*)m 

The  above  relationship  does  not  have  a  theoretical  justification  and  it 
will  be  considered  as  an  expediency*  At  any  rate,  (20)  reduces  to  i>  “  n , 
for  y  <  y*  and  if  =  V2*  for  y  >  y*,  and  to  ,  V*  =  {V* at  y*..  The 
arbitrary  exponent  5  tends  to  narrow  the  diffused  range  drastically. 

The  introduction  of  the  above  remedies,  however, does  not  conflict  with  our 
objectives,  because  the  correct  order  of  magnitude  of  eddy  transport  coef¬ 
ficients  are  maintained,  •  . 

Ill  all  computations,  the  following  values  of  .parameters  are  used; 
n  =  0.124,  K  =  0.4,  y*  =  2jV}]XJ  f  ,  aC=  l,  V  =  0,005  cm^/sec,  and  ; 

L  -  20  cm*.  Some  results  of  numerical  computation  will  be  shown  here  for 

P  =  0.008  and  C„At/X=  0.01.  In  Figure  3 »  the  turbulent  velocity  profiles 

r  ■  ■  ,  • 

at  X  s  Id  and  20  cm  are  compared  with  the  laminar  distribution.  The  effect 

of  eddy  kinematic  viscosity  is  to  flatten  the  velocity  profile  in  general, 

*3  The  particular  choice  of  values  for  V  and  L  is  immaterial  to  nur  results. 
In  particular,  the  appearance  of  V  could  have  been  eliminated  completely 
•  by  a  proper  transfopnatibn.  •. 
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The  lainAnsr  tenperature  proTile  is  essentially  linear  anu.  the  turbulent 
ones  seem  to  swing  around  the  linear  profile  vxith  small  deviations  as 
shown  in  Tigure  4.  In  Figure  5,  the  film  thicknesses  for  (a.)  Ilusselt’s 
lamitiar  case,  (b)  the  two-p'nase  laminar  case,  and  (c)  the  two-phase  tur¬ 
bulent  case  are  compared.  Because  of  the  existence  of  similarity,  (a) 
and  (b)  are  not  only  straight  but  displaced  by  a  constant  factor;  while 
(c)  deviates  from  (b)  and  increases  rapidly  with  x. 

Heat  Transfer  Results 

The  heat  transfer  coefficient  for  condendation  is  defined  as: 

(21)  h  »  k(St/dy)^j/ At 

which  is  obtained  from  the  heat  energy  balance,  hAt  =  k(bt  /Sy  In  , 
the  absence  of  similarity,  the  dependency  of  the  heat  transfer  coefficient 
on  the  plate  length  will  be  avoided  by  adopting  the  follovjing  parameters: 
Reynolds  number: 

,5 

(22)  *  hV /  -  (4/^)1  udy  at  x  ■  L 

'o 

and,  average  heat  transfer  coefficient: 

(23)  (  g  \y/  k  ■  (  g  (  1/At  L  )(  Ax 

*0 

For  Bussedt’s  laminar  case,  the  follovjing  relationship  between  Rex  and  h  „ 

(  -  1.U7  Rej, 


exists; 


(24) 
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In  Figure  6,  equation  (24)  »as  plotted  as  curve  (1)  and  the  results 

or  Husselfs  tnrhulent  oase-tata  Iron  rel.  W-v,ere  shown  by  curves  (2) 

and  (3),  Pr  =  analysis  are 

rll^  f'i'i  and  (7).  Curves  (4)  and  (5)  are  the 

su-nmarized  in  curves  (4),  l5J*  ana 

results  of  the  two-phase  lanlnar  case  for  P,  =  0.003  and  0.003  respectively. 
The  results  of  the  two-phase  turbulent  case  are  shown  by  curves  (6)  and  (7) 
for  the  sa.o  Prandtl  numbers.  As  in  the  previous  worh  [4],  the  heat  transfer 
results  are  not  sensitive  to  the  particular  value  of  cC  and  the  average  heat 
transfer  coefficient  is  decreased  by  not  wore  than  >4  by  taking  cC  =  0.1 

as  an  example » 


Discussions 

n  n  4-,^  o  non  which  covers  most  of  the  experi- 

Uithin  the  range  of  Re^  up  to  2,000,  whicn  covei 

nental  data  of  Hlsr.  and  Bonilla  t5L  curves  (2)  and  (3)  show  sowewhat  lower 

h  than  curve  (1),  but  the  difference  is  rather  nodest  (about  10?).  Sim- 

ulrly.  the  comparison  of  two  pairs  of  curves,  (4)  with  (S)  and  (5)  with  (7). 

indicates  that  the  turbulent  model  reduced  h^,,  roughly  by  the  same  fraotuon, 

i;e.  10?.  From  these  comparisons  one  can  conclude  that  the  turbulent  trans- 

ports  can  reduce  the  heat  transfer  coefficient  but  only  by  a  small  fraction. 

This  can  be  visualised  from  the  fact  ttot  the  increase  of  fil"  thickness  is 

counteracted  by  a  smaller  turbulent  velocity  so  as  to  give  approximately  the 

same  flow  rate  as  in  the  la.mi.iar  case.  Therefore,  any  decrease  in  the  heat 

transfer  coefficient  appears  to  be  directly  profKirtional  to  the  increase  in 

the  film  thickness,  because  h  ?  k/5.  However,  the  average  effect  fever  L)  is 


much  weaker. 


Nevertheless,  one  of  the  observations  vrorthy  of  attention  is  that 
curves  (4)  and  (5)  can  predict  considerably  lower  hay  corapared  to  curve 
(l).  This  is  certainly  attributable  to  the  improvenents  over  the  Nusselt’s 
model.  One  of  the  nodifications  notably  the  ”no.n-vanishing  interfacial 
shear  stress’’,  is  believed  to  be  mainly  respor slble  for  such  reduction  (see 
ref.  [9}  ).  All  in  all,  the  bulk  of  irregularly  scattered  data  still  lies 
well  below  all  the  curves,  therefore,  it  seems  logical  to  then  ask,  what 
would  be  the  effect  of  an  upward  vapor  velocity  in  liquid  metals  conden¬ 
sation.  This  was  motivated  by  the  fact  that  the  upvrard  interfacial  drag 
vre  have  considered  was  caused  by  a  quiescent  vapor  and  the  existance  of 
upward  vapor  flow  was  obsen,’'ed  by  Misra  and  Bonilla.  They  have  reported 
that,  in  some  of  the  jOv;_ pressure  film  condensation  runs,  the  vapor  flow 
was  large  enough  to  blow  off  the  film  from  the  condenser  plate  in  all 
directions.  Musselt  [5]  first  studied  the  effect  of  upward  vapor,  flows 
but  his  results  are  not  suited  for  oxir  discussion.  Therefore,  we  shall 
present  here  a  qualitative  analysis  on  the  effect  of  an  upvrard  vapor  flow 
based  on  Nusselt’s  model.  The  equation  of  motion  can  be  modified  as: 


(25)  du/dy  -  g  8  (1  -  y/8)/ 

where denotes  the  interfacial  shear  stress.  In  the  interest  of  main¬ 
taining  similarity,  let  us  set  Y- /  =Bg  S,  then  we  have: 

P 


(26) 


du/dy  -  g  8  (1  -  p)(l  -  y/8(l-j?))/  i 


where  o<^<2/3  is  assumed  (^<2/3  assures  r>®).  In  this  way,  the 
effect  of  upvrard  velocity  is  related  directly  to  the  position  of  zero  ve¬ 
locity  gradient  and  for  the  case  of  quiescent  vapor  of  Figure  3,  0,3 

(for  the  laminar  case).  A  detailed  investigation  of  condensation  using 


14 


(26)  is  rather  involved  because  the  character  of  solution  changes  as  the  inter¬ 


facial  shear  stress  increases,  Ilof-rever,  vdthin  the  range  cf 
velocity,  let  us  assume  that  the  boundary  condition  of  €  = 
still  valid,  i/ith  a  linear  tevajjerature  profile  and  constant 
sliip  betvreen  Re^  and  h^^,  siiailar  to  (24)  becomes: 

(27)  (  )  VA  -  1.47  (1  -  > 


small  upv/ard  vapor 
0  at  X  =  0. is  . 


At,  the 


relation 


ihe  above  equation  implies  that  h^^-  can  be  smaller  loy  a  factor  of 

1/3  ■  . 

(1-3p/2)  in  comparison  to  the  Nusselt’s  laminar  case  (24),  In  Figure  6, 
equation  (27)  vras  plotted  for  several  values  of  ^  ,  Curves  (C),  (9),  (lo), 
and  (11)  cover  most  of  the  data.  It  is  seen  that  the  upv/ard  vapor  flow  can 
cause  the  reduction  of  and,  consequently,  the  .scatter  of  data,  could  have 
been  caused  by  the  variations  in  vapor  velocity.  This  is  in  confimity  vlth  ths 


previous  findings  of  Colburn  and  Carpenter  [10],  As  one  expects,  hgy 
Siiould  increase  with  the  positive  (dovmward)  interfacial  shear  stress  [8]. 
However,  it  is  very  strange  to  note  that  the  data  of  Kisra  and  Bonilla  are 
not  consistent  vrith  the  trend  predicted  by  equation  (27),  In  fact,  they  have 


observed  that  the  data  seem  to  approach  equation  (24) ‘as  the  upward  vapor 
velocity  increases,  but  the  data  fall  vrell  below  (24)  when  the  upward  vapor 
flow  is  small.  They  pointed  out  that  these  facts  are  ’‘striking'’,  but  did  not 
Oi i er  an  explana tion.  It  must  be  noted  that  the  experimental  data  are  not 
very  accurate  due  to  the  unsurmountable  experimental  difficulties.  In  a  few 
instances,  the.  author  was  able  to  observe  much  higher  h  than  the  reported 
data  value  by  making  re-computations. 
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Conclusions 

In  an  attempt  to  explain  the  discrepancy  existing  betvjcen  the  theoretical 
predictions  and  the  experimental  data  for  liquid  metal  condensation,  the  two- 
phase  turbulent  condensation  vras  studied.  As  in  the  ‘lusselt’s  turbulent  case, 
the  inclusion  of  turbulent  transports  reduces  the  heat  transfer  coefficient. 

The  overall  effect  is  rather  modest  and  the  reduction  oi*  h^;^  does  not  exceed 
lOf-  over  the  laminar  model  within  the  experimental  range  of  In  experi¬ 

ments,  however,  the  presence  of  an  upward  vapor  flow  was  observed,  which  vras 
perhaps  inevitable  from  the  usual  experimental  set-up.  It  has  been  shown 
qualitatively  that  the  upward  vapor  flow  can  reduce  the  heat  transfer  coef¬ 
ficients.  Consequently,  the  scatter  of  data  could  have  been  caused  by  the 
variations  in  vapor  velocity.  In  order  to  confirri  the  assertion  proposed  , 
above,  it  is  necessary  to  carry  out  an  exact  fonaulation  of  the  problem  SO  as 
to  relate  r  to  the  vapor  velocity  directly  and  it  is  also  necessaiy  to 
measure  the  accurate  local  vapor  velocity  in  condensation* 
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Subscripts 

av  average  over  L 

i  interface 

s  saturation 

t  total  (molecular  -i-  eddy) 

w  v;all 

1  Deissler’s  region 

2  von  Karman’s  region 


Superscripts 

*  eddy  transport  coefficient 

/s  , variable  and  property  in  vapor  boundary  layer. 
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